MPP-2004-43 
May 2004 



Structure functions of 2d integrable asymptotically free 

models 



Janos Balog 

Research Institute for Particle and Nuclear Physics 
1525 Budapest 114, Pf. 49, Hungary 

Peter Weisz 

Max-Planck-Institut fur Physik 
Fohringer Ring 6, D-80805 Miinchen, Germany 

Abstract 

We investigate structure functions in the 2-dimensional (asymptotically free) non- 
linear O(ra) cr— models using the non-perturbative S-matrix bootstrap program. In 
particular the exact small (Bjorken) x behavior is exhibited; the structure is rather 
universal and is probably the same in a wide class of (integrable) asymptotically 
free models. Structurally similar universal formulae may also hold for the small 
x behavior of QCD in 4-dimensions. Structure functions in the special case of the 
n = 3 model are accurately computed over the whole x range for —q 2 /M 2 < 10 5 , and 
some moments are compared with results from renormalized perturbation theory. 
Some remarks concerning the structure functions in the \ jn approximation are also 
made. 



1 Introduction 



Structure functions describing scattering of electrons and neutrinos off nucleon tar- 
gets are well measured and give us insight into the structure of the nucleons jll2lol4j . 
At high — q 2 and intermediate Bjorken x the parton model and DGLAP equations 
|5l6j give a good description. At smaller x the DLGAP equations are considered 
to break down and BFKL-type |7j equations take over, here the structure function 
i*2 ~ x~ u (~ q ) with u(—q 2 ) > 0. A value of v(0) > would however (seem to) violate 
the Froissart bound. Recently saturation models, such as the so-called color glass 
model |8I9| predict the true asymptotic behavior to be ~ \v? x , with p = 1 or 2. 

The QCD literature on small x physics is vast and rather involved |10|llj . One 
certain aspect is that a description of the asymptotically small x region requires some 
crucial non-perturbative input. The most systematic non-perturbative methods for 
QCD, using the lattice regularization, are able to give non-perturbative information 
on the moments of the structure functions via the OPE |12|13j . however they are 
not applicable to yield information on the asymptotically small x behavior. 

In this paper we study structure functions in asymptotically free integrable mod- 
els in two dimensions. Despite the fact that there are no transverse directions, the 
structure functions have a rather rich and non-trivial behavior with many features 
reminiscent of the structure functions in QCD. Here we will concentrate on results 
obtained for O(n) sigma models. In particular we have derived the exact small 
x behavior in these models. The result has a rather universal model-independent 
structure, being independent of n and holding for a large class of operator probes. 
One is tempted to speculate a similar qualitative structure to hold in QCD. 

2 Sigma model 2— point functions 

The O(n) u-model in 2d formally described by the Lagrangian 



is perturbatively asymptotically free for n > 3. A special property is that these 
models have an infinite number of local Jl] and non-local JH] classical conservation 
laws which survive quantization. At the quantum level they imply absence of particle 
production. Assuming the spectrum to consist of one stable 0(n)-vector multiplet 
of mass M, the S-matrix has been proposed long ago by the Zamolodchikovs |16j . 
Form factors of local operators can be computed using general principles |17| . The 
S-matrix bootstrap program for the construction of correlation functions involves 
summing the contributions over all intermediate states |18] . The possible equivalence 
of this construction to the continuum limit of the lattice regularized theory has 





1 



been investigated in ref. [T^]. In papers of one of the present authors (J.B) and 
M. Niedermaier |2U| 2-point functions of various operators were computed, including 
those of the O(n) current J^f and spin-field operators <3? a . We give their definitions 
here because they will be needed later to state our result on the small x behavior: 



(0\TJf(x)J e J(y)\0) 

A 2 

_i (v*sF - 8*5*1*) I 9_ e -i*(*-„) _ q 2 v ^ h{ _ q 2 _ - e) ( (2 2) 



valid up to contact terms and 



A 2 

(0\T<5> a (x)$ b (y)\0) = -i&*A* J e -*«(*-v) ^(.^ _ ie) , (2 . 3 ) 
where the normalization factor A n is chosen (for later convenience): 

A 3 = -^=, A« = l n>3. (2.4) 



To complete the definitions we must supplement the operator normalizations. The 
currents have a normalization fixed by requiring their spatial integrals to yield the 
correct charges (cf Eq. (|4.11|) ). and our field normalization is fixed by requiring 

(0|$ a (0)|&,6>) =S ab (2.5) 

for one particle states with momentum p\ = Msinh#, with state normalization 

(a,0'\b,6) = 4ir5 ab 5(9' -8). (2.6) 

The studies of these 2-point functions (in the case n = 3) [2U] presently constitute 
the best evidence for the existence of a non-perturbative construction of a model 
with asymptotic freedom. 

3 Sigma model structure functions 

The (T-model analogue of the central object in deep inelastic scattering is 

WfiMfaq) = ^J d 2 xe^(a,p\[jf(x),J e J(0)]\b,p). (3.1) 

We define the usual DIS kinematic variables 

u=pq/M, x = -q 2 /{2Mv) , W 2 = (p + qf . (3.2) 
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In the relevant kinematic domain q 2 < 0, W > M i.e. < x < 1 we have 

^ ;cde/ (p,g) = ^^<a,p|J^(0)|r)(r|^(0)|6,p)^ 2 )( j 9 + tZ - i^.), (3.3) 

r 

where the sum is over the complete set of r-particle ("in" or "out") states. Using 
Lorentz and O(n) invariance we have 

W^(p, q) = Uu ~ £ vntf, x)R?^ , (3.4) 

v q J 1=0 

with projectors Ri corresponding to the 3 invariant i-channel "isospins" given by 

R ab;cdef = 1 ^6 sgce $df _ fief gdey ,35) 

n 

Rf cdef = ^(S ac S be -S bc 5 ae )S df -(c^d)-(e^ f) + (c~d,e~ /), (3.6) 

R ab;cdef = _ 1 j ^ § be + flogae^ _ ( c <_> ^ _ ( e <_> y) + ( c <_> d ^ & <_> J) j 

+ 2 £a6 ^c/^deN (3.7) 

n 

Note in 2 dimensions there is only one structure function for each isospin channel, 
since there is only one (conserved symmetric) tensor involving two momenta; one 
has e.g. ( for p 2 = M 2 ) 

/ Muq^ ( Mvqy\ M 2 {u 2 - q 2 ) / 2 x 

" -^2" J [P» ~ ) = ^2)2 " ^9 J • ( 3 - 8 ) 

Although in QCD structure functions for operators not associated with physical 
sources have so far not been studied, we also introduce, for reasons which will soon 
become apparent, the field structure functions through 

Z ab -> cd (p,q) = -ir q 2 J2(a, P mO)\r)(r\3> d (0)\b,p)S {2) (p + q-Pr) (3.9) 

r 

= A 2 n j2MQ 2 ,x)P^ cd , (3.10) 

1=0 

with t-channel projectors Pi for the vector representation given by 

pab;cd = ]_ 5 ab 5 cd^ (3 n ) 



n : 2 



}_^ S ac S bd _ S bc 5 ady ( 312 ) 
pafe;cd = l/ (5 ac (5 6d + (5 &c (5 a^ _ I^cd (3.13) 
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The current operators connect only states with an even number of particles to the 
vacuum and the field operators only states with and odd number: 



wi(q 2 ,x) = wl r \q 2 ,x) , w t (q 2 ,x) = ^ w { [\q 2 ,x). (3.14) 

r odd r even 



The wonderful feature of integrable models is that one knows many explicit 
properties concerning the form factors appearing in the expressions. These are en- 
capsulated in the Smirnov axioms , and using these one can derive the previously 
announced exact asymptotic small x behavior of the structure functions: 

Wliq2 > x) ~ xll 2 x(n- T 2)i M - q2) > (3-15) 

™ i{q2 > x) ~ xL 2 x{T-^ M - q2) ' (3 - i6) 

where the functions A s are the "Adler functions" defined from the vacuum 2-point 
functions (|2T2j) . l|23)l by 

A s {z) = -z 2 ^-I s {z), « = 0,1. (3.17) 

The constants appearing in (|3.15|) . ()3.16|) are characteristic of the O(n) representa- 
tions (vector V and anti-symmetric tensor T) of the corresponding operators: 

V = 2(n - 1) , V 1 = V 2 = n-2, (3.18) 
T = 4(n-2), Tx = n-2, T 2 = 4 - n . (3.19) 

The explicit proof of these relations will be presented in a forthcoming more detailed 
publication In this letter we would like to concentrate on general features. 

Firstly we note that the structure of the asymptotic small x behavior is independent 
of the operator, independent of n, and independent of the channel. Further since 
the results were obtained from rather general principles, we think that they are valid 
for a large class of integrable asymptotically free models. Indeed we have checked 
that the same behavior holds in the leading orders of the 1/n expansion in the 
Gross-Neveu model |22| . 

In the derivation of (|3.15[) . (|3,16|) it appears that the 1/ In 2 x behavior is related 
to the high energy behavior of the scattering amplitudes; this is similar to the 
association of the proposed ln p x in QCD with the asymptotic behavior of total cross 
sections (related to the forward scattering amplitude through the optical theorem). 

The question is what can one learn for QCD; is the small x behavior there also 
given by a structural formula factorizing a part characteristic to the target and a 
part described by the vacuum 2-point function? Unfortunately so far we have not 
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Figure 1: F2(q 
with error bars), 



,x)/D(-q 2 ) for x = 5 • 10 
for various values of — q 



(circles with error bars) and x = 8 • 10 a (squares 



2 (in GeV 2 ). Symbols without error bars are values of 



xwo(q ,x)/Ai(—q ) in the 0(3) <r-model at the same values of x (here the abscissa denotes values 
of -q 2 /M 2 ). 



succeeded to derive such a general result. As a first guess we have looked at the 
ratio of the structure function F<i to the (electro-magnetic) Adler function D in 
QCD, with HERA data at some of the lowest x values published so far |24| . and the 
Adler function taken from ref. [23]. The result is presented in Fig. 1; there is no sign 
that the ratio is becoming independent of — q 2 as x decreases. However from such 
comparisons one should be cautious to draw conclusions concerning the asymptotic 
behavior, because at these values of x the situation is qualitatively similar (only the 
slope is different) for the ratio xwo(q 2 , x)/A\{ — q 2 ) in the 0(3) <7-model, which is 
also plotted in the same figure. In this model one would have to go to much smaller 
values of x. The question remains for QCD, at which value of x (for a given — q 2 
range) does the true asymptotic behavior set in? 

The QCD structure functions in the range of small x between 10~ 5 and 1CP 2 
are fitted quite well with a "Lipatov form" A(— q 2 )x~ v (~ 9 \ As an exercise we have 
made least-square fits of xwo(q 2 , x) (for 0(3)) with such a form and observed that 
in the regime 1(T 5 < x < 1(T 3 such fits with v(-q 2 ) = -0.192, -0.167, -0.151 for 
—q 2 /M 2 = 1, 10, 100 respectively, describe the data better than simplest fits of the 
form C(—q 2 )/ In 2 x incorporating the known asymptotic small x behavior. 
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3.1 1/n expansion 



The data obtained for 0(3) used above and in subsequent sections require a consid- 
erable amount of computation. For this reason we include here a short discussion of 
structure functions in leading order of the 1/n expansion, where many qualitative 
results are rather similar to n = 3 but where relatively simple analytic formulae are 
available. The easiest case is the leading order of the spin structure functions in the 
isospin 1,2 channels because they are given by the imaginary part of the propagator 
of the "auxiliary field": 

i _ rfi /If 2 sinh B 

w^^x) = -9(W - 2M)4vr { _^ + m2? q2 - ^ + 0(l/n 2 ) , 

for W 2 = AM 2 cosh 2 °-, 1 = 1,2. (3.20) 

Apart from suggesting that the limits n — > oo and x —* commute, this simple 
function already illustrates many rather general features of the structure functions 
in this model. Firstly that the onset of the limit x — > is not uniform in — q 2 . In 
the Bj-limit — q 2 — ► oo , x fixed we have 

~ / •> 1 2tt(1 — x) , 

w ^\x) ; 27 2 I L^ l = ^ 2 - 3 -21 

nxm (— q 2 /M 2 ) 

Secondly the limit x — > for fixed q 2 is approached extremely slowly e.g. for 
-q 2 /M 2 = 1 , wi/(wi) asympt = 0.93 for x = 10" 5 , while for -q 2 /M 2 = 100 we have 
e.g. i^z /(i'i)asympt = 0.49,0.59 for x = 10 -5 , 10~ 7 respectively. 

We remark that the leading order (in 1/n) isospin (field) structure function 
involves another (box) diagram which is also easily evaluated. One can show that 
the small x behavior is as predicted by the general formula (|3.16|) . The leading 1/n 
diagram contributing to the current structure functions is just a 1-particle exchange 
and thus only contributes a term oc 5(1 — x). We have computed the next-to-leading 
orders for the I = 1,2 channels and again confirm the predicted behavior. 



3.2 The case n = 3 

So far we have concentrated on the small x region; in the following we extend our 
description of the structure functions to the whole range of x. We will do this for the 
case n = 3 which is rather special for various reasons. For our studies in the S-matrix 
bootstrap approach it is important that it is the model for which the r-particle form 
factors can most easily be obtained explicitly. Moreover the spin and current 2-point 
functions exhibit in this case very similar features and there are miraculous scaling 
relations [20] which relate them . There are well defined recursive procedures for 

1 see also the OPE in sect. 4.2 
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Figure 2: Approximations to a:(uio + wo) as functions of < x < 1 for —q 2 /M 2 — 100. Curves 
correspond to sums up to and including 2, 3, 4, 5, 6-particle intermediate states. The last 3 curves 
are indistinguishable on this scale. 

computing the form factors, the only practical limitation being that they become 
very involved. So far the record we have achieved is the 7-particle form factor of the 
spin field [25]; already its algebraic expression in MAPLE involves many megabytes 
of storage. Fortunately for the structure functions we only require sums over bilinear 
factors of the form factors which are computationally more manageable. 

Just as for the 2-point functions we find that for a fixed — q 2 only states 
with a limited number of particles contribute significantly. To appreciate this better 
we consider the sum of the field and current structure functions, which is a rather 
peculiar thing to do in general, but which is in fact rather natural in the special 
case n = 3. Figs. 2 and 3 illustrate how the structure function x(wo + wo) is built 
up from states with increasing particle number for the cases —q 2 /M 2 = 10 2 and 
—q 2 /M 2 = 10 4 respectively. We see that the higher states contribute very little 
and that we obtain nearly exact values for the structure functions for all values 
of —q 2 /M 2 < 10 5 by including only intermediate states with < 5 particles for the 
current and < 6 particles for the spin field. 

In Fig. 4 we plot xwo(q 2 ,x) as a function of log 10 (— q 2 /M 2 ), for a selection of 
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Figure 3: Approximations to x(wo + wo) functions of < x < 0.95 for —q 2 /M 2 — 10 4 . Curves 
correspond to sums up to and including 2,3,4,5,6-particle states. 



x-values 2 . The function increases as — q 2 increases for all values of x in this range. 
3.3 Threshold behavior 

Note that in Fig. 3 we have cut off the plot at x = 0.95. This is because near x = 1 
the function develops a big bump with a peak ~ 70 which, if included in the same 
plot, would obscure the features we wanted to show there. The behavior of the 
c-model structure functions near x = 1 is indeed rather involved. For a fixed —q 2 
the contribution to the structure function from the r-particle state vanishes for 
x greater than some threshold value 

x r {-q 2 ) = [l - (r 2 - l)M 2 /q 2 ] ~* . (3.22) 

The big bump in x(wo + wo) referred to above is at this value of —q 2 /M 2 = 10 4 
practically entirely due to the 2-particle contribution. For this contribution: 



wf ] ~ E l {q 2 )Jx 2 (-q 2 )-x, x x 2 , -q 2 fixed, (3.23) 



2 For this model we prefer to show this rather than the typical HERA plot where one adds 
— log 10 (a:) to separate the x-values, because the latter would obscure the — q variation which is 
rather small compared to the variation of — log 10 (a;). 
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log_10(-q A 2/M A 2) 
Figure 4: xw (q 2 ,x) for various values of x = 1CP 1 / 5 . 

wf ] ~ Fi(x)/ In 2 (-q 2 /M 2 ), -q 2 -> oo , x fixed, (3.24) 

where Ei,Fi are some (known) functions. The bump arises because Fi is quite 
singular near threshold, Fi ~ [(1 — x) ln 2 (l — a;)] -1 . The analytic behavior as x — ► X2 
sets in only extremely close to threshold e.g for —q 2 /M 2 = 10 4 the position of 
the peak of the bump is at x = 0.99954 whereas the function vanishes at X2 = 
0.99970. At —q 2 /M 2 = 10 4 the 3-particle contribution also has a bump but it is less 
pronounced; (peak value ~ 2.5 at x ~ 0.9953). We conjecture that the threshold 
behavior of vJlp in the 0(3) model is (x r — x)^' 2 ^/ 2 . 

One can also study the threshold behavior in the n — ► oo limit. Here we find 
that in the leading order wi ~ \/x2{— q 2 ) — x for I = 1, 2 but u>o ~ l/^/x^— q 2 ) — x. 
We caution however that the limits n — > oo and x — > threshold may not commute. 

The characteristic enhanced near-threshold features in this model probably 
have no counterpart in QCD. For QCD the behavior of the structure functions 
in the x ~ 1 region is surely also complicated. But there the effects might be 
quite suppressed for large — q 2 since, if we model the behavior in this region by the 
contribution of resonances, their electromagnetic form factors are thought to fall 
very quickly (as powers) in — q 2 similarly to that of the proton. 
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4 Partons, OPE and moments 



4.1 Parton model 

In the O(n) a-models there does not seem to be a simple parton picture. This is 
even so for the case n = 3 where the model is equivalent to the CP 1 model. For 
although this model is formulated in terms of a complex doublet of fields which 
are analogous to quarks in that they are confined, it seems that they do not play 
a role more similar to partons than the elementary bare spin fields in the original 
formulation 3 . The question is related to that of understanding what are (if any) 
the "ultra-particles" in the sense of Buchholz and Verch |26j . or to the associated 
question as to whether the c-models have an underlying conformal field theory. 
Although an intuitive parton description with suggestive DGLAP equations 

d f 1 dy 

Q 2 Tr^wi(q 2 ,x) = — pi{x/y,q 2 )wi(q 2 ,y) , (4.1) 
oq z J x y 

(where pi(z,q 2 ) would be the corresponding splitting functions) is still missing in 
these models, we still have the machinery of the operator product expansion (OPE) 
which we apply in the following. 

4.2 Moments 

A class of interesting quantities are the moments of the structure functions: 

E M kN^) , (4-2) 

r=odd 

E iO 2 ), (4.3) 

r=even 

where M^{q 2 ) denote the contributions from r-particle states 

Ml 7 i(q 2 )= fda^-^V.*). M) 
J o 

As in QCD the moments satisfy renormalization group equations from which one 
can determine their leading behavior as — q 2 — > oo. These are derived by considering 
the OPE for the current-current and field-field products and explicitly treating the 
terms involving the operators of highest (zero) twist 4 . The general analysis is 

3 Perhaps the peculiar threshold behavior discussed in sect. 3.3 is explained by the fact that (as 
opposed to QCD) with some probability the O(n) particle can consist of a single point-like parton 
that carries the same quantum numbers. 

4 "Twist" in this model is the naive dimension minus the "spin" (^number of uncontracted 
Lorentz indices). 



M hN (q 2 ) = [ dxx N ~ 1 w l (q 2 ,x) = 
Jo 

M kN (q 2 ) = [ dxx N - 1 wi{q 2 ,x) = 
Jo 



10 



rather involved because the classification of lowest twist operators turns out to be 
more complicated than in QCD because the elementary field is dimensionless 5 . Our 
analysis extends that initiated e.g. in refs.|15|. [2*7] . Here we just quote some results 
and defer the derivations to ref. |21| . 

For the current (N even) moments in the isospin channel we have 

M Mq 2 )=W ., N -^- \l + ^—X(q 2 )+0(\ 2 )} , N>2, (4.5) 
2{n —1)1 n — 2 V/ J 

where X(q 2 ) is an effective running coupling function defined through 



' +^_lnA(g 2 )=ln^-^, (4.6) 



X(q 2 ) n-2 ■ A m 

and the Wo-n are renormalization group invariant, non-perturbative constants, cor- 
responding to the matrix elements of spin N operators. In the N = 2 case this is the 
energy-momentum tensor operator T^ v for which we know the constant explicitly 

(a,e\T ++ (0)\b,9) = - A W , 2 M\ 2e 5 a \ W -2 = 2 , (4.7) 

where the index + means "(0 — l)/2". In particular the "momentum sum rule" 
follows: 

n — 2 

M 0;2 (-oo) = . (4.8) 

n — 1 

Note that all the isospin moments tend to constants as — q 2 — > oo. As a consequence 
these current structure functions in the O(n) models obey Bjorken scaling, and 
Fig. 4 indicates that the resulting limiting scaling functions are non-trivial. This is 
a special property of these models and we conjecture that this is due to the existence 
of an infinite set of local conserved quantities |14j . 

In the isospin I = 1 channel for odd moments N > 3 we can only say that 

M 1;JV (g 2 ) = W 1;N X(q 2 )^ + . . . , N > 3 , (4.9) 

but in the special case TV = 1 we have 

M 1 , l (q 2 )= l - |l-_l_A(g 2 )+0(A 2 )| , (4.10) 

where the constant is known through the current normalization 

(a,0\jf{Q)\b,6) = -2iMe e Pf ;cd . (4.11) 

5 cf. in QCD the quark field carries dimension 3/2 
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From this follows the analogy to the Adler sum rule in QCD: 



2 

For the spin field isospin moments we have 



Mi ; i(-oo) = ^. (4.12) 



M o;iv(9 2 ) = { w { z: 2 J . „; _ \.L ~ . (4-13) 



0:.V 



{1 + A(g 2 ) + 0(A 2 )} 



where the non-perturbative constants Wo-N are the same as for the current, and 
where C n is the non-perturbative constant appearing in the short distance expansion 

(0|$ a (y)$ 6 (0)|0> ~ C n 5 ab (-lnAf|y|)^ , (4.14) 

which is only known for the cases n = 3 and n = oo (Coo = l/(2vr)). We see that 
only for the case n = 3 do the moments of the field I = structure function have 
the same leading asymptotic behavior as those of the current. 

For the isospin I = 1 field (odd) moments we find to leading orders PT 

~ 2 fM , 2 (q 2 ), N=l; 

Ml-NW ) = { ~ 2n-5 (4.15) 

I W^nX— + . . . , N > 3, 

where there is in general no obvious relation between the W\-n and the constants 
occurring in l|4.9jl . except for n = 3 where they are equal (Wi-n = W\-n , n = 3). 

(r) (r) 

In Figs. 5 and 6 we plot the separate r-particle contributions Mq. 2 and M^.-[ 
respectively (for n = 3). They are typically bell-shaped (except for r = 1) and 
perhaps obey scaling relations similar to those of the spectral functions examined in 
ref. [2U]. The figures show how they build up the sum of moments Mq-2 + Mq-2 and 
Mi ; i + Mi-i. Using the exact ratio of the mass to the A-parameter 

M ( 8 / e )i/(n-2) 

A=T- r [l + l/(n-2)]' ( - 16) 

obtained by Hasenfratz, Maggiore and Niedermayer [2Hj) we also exhibit the per- 
turbative results up to and including terms of order X(q 2 ). The agreement of the 
summed terms and PT is impressive for —q 2 /M 2 ~ 10 5 . For values of —q 2 /M 2 >~ 
10 6 contributions from states with > 7 particles must be taken into account. Note 
we have also included the contribution of the one particle states in the sums; these 
tend to improve the agreement at lower values of — q 2 and fall asymptotically as 
Af/.JJ ~ miTT 4 /[Aln 2 (-q 2 /M 2 )} with m = l,mi = 1/2, m 2 = -1/2. 
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Figure 5: Contributions Mg r 2 for n — 3 from r = 1, . . . , 6-particle states. The upper full line is 
their sum. The dashed line is the perturbative expansion of Mo-2 + Mo-2 = 1 + A up to and including 
terms of order \{q 2 )- 



4.3 Conclusions 

Many qualitative field-theoretic features first observed in non-perturbative studies 
of integrable models in 2d, have in the past found their counterparts in realistic 
models in 4 dimensions. Although fascinating in their own right, we hope that the 
investigations of structure functions of 2d O(n) u-models presented in this paper 
will play a similar role. Similar methods are applicable to many other physical sit- 
uations e.g. generalized structure functions, exclusive electro-production processes 
and rapidity gaps. 
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Figure 6: As for Fig. 5 but for the moment I = 1, N = 1; here the PT result is 1 + 0(A 2 ). 
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